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We harnesses the freedom in the celebrated gauge transformation approach to generate dark 
solitons of coupled nonlinear Schrodinger (NLS) type equations. The new approach which is purely 
algebraic could prove to be very useful, particularly in the construction of vector dark solitons in the 
fields of nonlinear optics, plasma physics and Bose-Einstein condensates. We have employed this 
algebraic method to coupled Gross- Pitaevskii (GP) and NLS equations and obtained dark solitons. 
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I. INTRODUCTION 

Eventhough the discovery of solitons in the numerical simulation of Korteweg de Vries (KdV) equation by Zabusky 
and Kruskal [I] goes back by almost fifty years and several nonlinear integrable partial differential equations (pdes) 
like modified KdV (mKdv) [2], Sine Gordon(SG) [3], nonlinear Schrodinger (NLS)[T equations have been identified, 
construction of integrable dynamical systems and the associated soliton solution continues to be a challenging problem 
even today [5]. In the case of NLS equation, soliton type pulse propagation in the anomalous and the normal dispersion 
regimes are called bright and dark solitons, respectively[6-9]. The identification of bright solitons in NLS equation 
arising by virtue of the subtle balance between the Kerr nonlinearity and anomalous Group Velocity Dispersion (GVD) 
has made quite a turnaround in the fields of optical fibre communications and responsible for several experimental 
demonstrations mmmm- In addition to bright solitons, it was also reported that optical solitons can also be 
generated if GVD is negative (i.e., normal dispersion regime) leading to the possibility of different type of solitons called 
dark solitons 0E]. The use of optical solitons opens new prospects for non-interference transmission since solitons 
are highly stable with respect to perturbations caused by fibre nonuniformities and to external interference |7, 8, 12]. 
During the past four decades or so, it should be mentioned that the NLS equation has also been used to analyse the 
dark solitons. Recent studies on the dark solitons have revealed very interesting properties which may allow their 
stable transmission with much less spacing between solitons when compared with bright solitons. Also, the interaction 
effects between two dark solitons are less than the bright solitons in the presence of fibre loss. The interaction forces 
between two dark solitons are always repulsive, unlike the case of bright solitons where the interaction forces change 
according to their relative phase. The self-induced Raman effect is found to be more destructive in the case of dark 
solitons. Use of dark solitons for high-speed communication systems will remain an interesting subject for future 
research 00II31IT1]. It is worth pointing out at this juncture that several analytical techniques like Hirota bilinear 
method Hi cog, Darboux transformation method PH ng, inverse scattering method m have been employed to 
generate bright and dark solitons. 

An iterative approach by gauge transforming the eigenfunctions of the associated linear eigenvalue problem was 
developed by Ling Lie Chau et.al., m to generate soliton solutions of integrable nonlinear pdes. Eventhough this 
method which is quite similar to Darboux transformation approach is purely algebraic in nature as it enables one 
to construct soliton solutions from a trivial/nontrivial seed solution, the above approach has been employed only for 
generating bright solitons for zero seed solution and kink solitons for non-zero (constant) seed solution of nonlinear 
partial differential equations. In this paper, by choosing a plane wave as the seed solution, we suitably harness the 
gauge transformation approach to generate dark solitons of coupled NLS type equations. 

This paper is organized as follows: In section 2, we derive the algorithm to generate dark soliton solutions 
by playing around the freedom in the Gauge transformation method. In section 3, as an application of the above 
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algorithm, we experiment with coupled Gross Pitaevskii (GP) and coupled NLS equation and obtain the corresponding 
dark soliton solutions.We then discuss the collisional dynamics of dark solitons. 


II. MODIFIED GAUGE TRANSFORMATION APPROACH 

We know that any vector integrable (1+1) dimensional nonlinear Schrodinger equation (in general) can be described 
by the celebrated AKNS type linear eigenvalue problem of the following form 


(la) 

(lb) 

where U and V represent Lax pair governed by 3x3 matrices and = (</>i, 02,03) T is the associated eigenfunction 
of the linear eigenvalue problem so that the compatibility condition 


Ut-V x + [U,V]= 0 (2) 

generates the desirable integrable vector NLS equation. It should be mentioned that beginning with non-zero plane 
wave as the seed solution is crucial for generating dark solitons. Dark solitons as is well known, are dips in density 
profile IV’I 2 - And inorder to get such dips, one needs to have a non-zero background density to begin with which is 
ensured by the choice of non-zero plane wave as the seed solution. We now feed the nonzero plane wave solution into 
the Lax pair matrices V and V and obtain the associated linear system. 

$0 = U°$°, 

where 

^ == ^ | nonzero planewave seed • 

We now gauge transform the eigenfunction <f>° such that <f> = g&° where g is a gauge function represented by a 
3x3 matrix while is an iterated eigenfunction. The Lax representation in terms of the iterated eigenfunction is 
given by 

+ u 1 $ = 0, (5a) 

= 0} (5b) 

so that 


(3a) 

(3b) 


U^gU^g 1 I g x g 1 
Vi=gV°g~ l +g t g~ 1 


Now, the gauge function g(x, £), can be chosen in such a way that it represents the solution of the associated Riemann 
problem and it is meromorphic in the complex £ plane as 


g(x,t-,0 


1 + 


• J. 


(6) 


In the above equation, £ is the eigenvalue parameter while (i and (j represent arbitrary complex parameters and P 
is a 3 x 3 projection matrix. It should be noted that for scalar NLS type equations, J = a 3 (Pauli’s spin matrix) 
while for vector NLS type equation 


J = 


1 0 0 
0-10 
0 0-1 


(7) 



3 


It should be mentioned that, for scalar NLS equations, the gauge function g(x,t), U, V and the projection matrix P 
are (2x 2) in nature and <fi is a (2 x 1) column vector. The inverse of matrix g is given by 


9 = J ■ 


I -7 


( 8 ) 


Now, choosing £1 = £*, the projection matrix can be determined by solving the following set of partial differential 
equations 

Px = (1 - P)JU^(Ci)JP - PJU {0 \(i)J(l - P), (9a) 

Pt = (1 - - PJVM (Ci) J(1 - P), (9b) 

It should be mentioned that equations. © are identical to the one proposed by Ling Lie Chau et.al..[20] in the gauge 
transformation approach with a difference that P, J, and are (2x2) matrices for scalar nonlinear pdes while 
one has to work with (3 x 3) matrices for vector nonlinear pdes. One can then relate the projection matrix P(x, t) 
with the vacuum eigenfunction by 


where 


and 


P = J • P • J, 

p- M(1) 

Trace [MP)] ’ 

ilf(i) = $ (0) (^,i,Ci) -m (1) -$ (0) (a;,i,Ci) _1 . 


( 10 ) 

( 11 ) 


( 12 ) 


— ^(^? Cl)\nonzero planewave seed • 

In the above equation, /b ( 1 * is a 3x3 arbitrary matrix of the form 




# = 


e 2Sl y/2 e W e 2i (» + h) £ ( 1 ) e 2 i(Xi+? 2 ) \ 


r *(l)„-2i( Xl +€ 1 ) p—2<5i 


-*( 1 )p-2i(xi+?2) 


/V 2 


0 


0 

s- 251 /^ 


(13) 


(14) 


where are arbitrary functions of (x,t) and their choice is governed by the dispersion relation of the 

associated nonlinear pdes while and are the coupling parameters. It should be noted that the structure of 
m W matrix determines the nature (bright or dark) of soliton solutions of vector NLS type equations. 

Hence, one can write down the dark soliton solution as, 


d 1} = d 0) - 2i (Ci - Ci)A 2, (is) 

V>2 1} = V’2° ) - 2i(Cl - Cl)Pl3* (16) 

where ipf* ^(i=l,2) represents the seed solution while (i=l,2) denotes the iterated dark soliton solutions of the 
corresponding coupled NLS type equation and 

p Mk 

Mi i + M22 + M33 


p Mk 

Mh + M 2 >2 + M\ 3 

The modified gauge transformation approach can be extended to generate multidark soliton solutions. For example, 
the general form of ‘N’th dark soliton solution can be written as 


ki N) = ki N 1} - 2i(6v - Civ)-^, (17a) 


(17b) 
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where Pi 2 and Pi 3 are given by 


and 


6 AT -1 _ 
r 12 — 


?N-l 

13 


[Mg'~ 1} ((t + 7 M 1 ( f“ 1) )M 1 ( f } + 7 

M 2 ( f ) 7 7r 2 + m [% _1) A^f } )) + 

((r + + 7 (M 1 ( f- 1) M^ ) 

)M- 


M\ 


x n 

(JV-l) ,,(iV) 


^ VJ ))777 + ((r + 7 M 1 ( f- 1) 


W 

12 


r(N-i) M ( N ) 


7(-^12 1Ki 22 T J «13 

(r + M 2 ( f _1) 7 *)/r 2 ], 




= - [M 1 f- 1) ((r + 7 M 1 ( f- 1 ) )M 1 ( f ) + 7 (M 1 ( f- 1) 


72 

(TV-1) 

23 


M 2 ( f } + M 1 ( f- 1 ) Mif ) )) 7 *)/r 2 + M, 

((r + 7 M 1 ( f- 1) )M 1 ( f + 7 (C _1) ^ 


ilP 


71 

(TV-1) ..(TV) 


M7 ; ))777 + ((r + 7 M 1 ( f- 1) 


)M- 


(JV) 

13 


■M- 




(N-l) n AN) 


13 


V77 


T = _1) + m£~ 1] + M^~ 1} 


22 


73 


R = PR- 1 + P&- 1 + Pg** 


33 ^ 


7 = 

_ 


Cl - Cl 
C 2 - Ci ’ 
Cl - Cl 
C 2 -CV 


with 


Pfr 1 = l M 2 ( f- l) ((r + 7 M 1 ( f- 1) )M 1 7 ) +7(M 1 ( f“ 1) M 2 7 ) 

+ ^<>)) 777 + + 7M7” 1 )) 

+ m7- 1 ) M 3 7 ) ))77t 2 M 1 7 ) +7(M 1 7“ 1) 

+((r + 7 M<f-«)M<f) + 7 (eH* 1 + 

P^- 1 = M 1 ( ^“ 1) (7M 1 ( f ) M^f“ 1) +M 2 ( f ) (r + 7M^“ 1) ) 

+ 7 M^- 1 ) m7 ) )77t 2 + M 3 ( 2 Ar - 1 ) ( 7 Mi 7 ) M(f- 1) 

+ (r + 7 M 2 ( f “ 1} )M 2 ( f } + 7 M^“ 1) M^ ) ) 7 7 r 2 
+ ( 7 M 1 ( 2 Ar) M 2 ( f- 1) + (r + 7 M<7* ) )J< r) + 7 M^- 1) 
M 3 ( ^ ) )(t + M 2 ( ^“ 1} 7*))/t 2 

77 1 = M 1 ( 7 1 t7M7 ) M7- 1 77M7 ) M(7 1) +Mt Ar) 

(r + 7 Mif- 1 ) )) 7 7 r 2 + M 2 ( f- 1 ) ( 7 M 1 7 ) M(f- 1) 

+ 7 + M$\t + 7Mif- 1) ))7*/r 2 

+ ( 7 M 1 ( 3 Ar) M 3 ( f- 1) + 7 + (r 

+ 7 M 3 ( 3 Ar - 1 ) )MW)(r + Mjf-VVr 2 , 


Aijf 0 = e-^72; M 1 7 ) = e _ ^ A,+Xl ej JV) ; s e-^ef; 

.u 2| v) = r's.v • M^ ] =e eN /V 2; = 0; 

M^f 0 = e i?JV+Xl e 2 (iV) ; M^ } = 0; A/V° = e^/V2, 
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To generate dark solitons for the scalar NLS type equations, one can feed the same nonzero plane wave solution 
and follow the above procedure except that one has to choose the m W matrix of the following form 


m»> = ( "*■ }{"' 

n\ l/mi 


(18) 


k + iX 

rrii = —:—-— and n\ = 2 , 


(19) 


Then, equation. (15) yields the dark soliton solution for the corresponding integrable scalar NLS type equation. 


III. APPLICATIONS 


(i) Coupled GP (CGP) equation: 

The coupled GP equation representing a binary BEC comprising of the hyperfine states of the same atomic species 
say, Rubidium ( 87 Rb) in a transient harmonic trap can be written down (in dimensionless form) as 


, dy>j d^i 

dt dx 2 


2 

(t)'^2 b jk\i’k\ 2 +v{x,t) 
k =1 


V’i 


( 20 ) 


In the above equation, v(x,t) = A (t) 2 x 2 represents the transient harmonic trap and a(t) represents the temporal 
scattering length (a(t) should be negative for defocussing (attractive) and positive for focussing (repulsive) cases 
respectively) and j = 1,2 describes the order parameter of the condensates. The above integrable coupled GP 
equation has also been investigated [22] [23] and the dynamics of the vector BECs has been explored by constructing 
bright and dark solitons. 


The above equation. (20) admits the following eigenvalue problem, 


where, 0 = (</>i, 02 , foV and 


&X — 

U(j) 



4>t = 

Vcb 


(21) 


Qi 

Q2 \ 


= [ -Q1 - 


0 

(22) 

\-Q*2 

0 

-m J 


/Ui 

Vl2 

Vvi \ 


V = l v 21 

v 22 

^23 , 

(23) 

\ V 31 

V 32 

v 3 / 
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where, 


with 


Vn 

V 12 

V13 

V 21 

V 22 

V 23 

V 31 

V 32 

V 33 


-i((t ) 2 + ir(t)x((t ) + ~7(t) A (t)QiQi + l -l(t)A(t)Q 2 Q*2 
[r (t)x - C(t)]Qi + ^ Qix 

[r(i)x - C(t)]Q 2 + I Q 2x 

-[T(t)x-at)]Q * 1 + ±Q* lx 

i((t ) 2 - iT(t)x((t) - % ~7(t)A(t)QiQl 

-\Q 2 Ql 

-[T(t)x-at)]Q*2 + \Q*2 x 

-\QlQl 

i((t ) 2 - iT(t)x((t) - l -l{t)A(t)Q 2 Q * 2 


Qi- ^=M^y { ~ T{t)x2 \ 
\fW) 

Mx,t)e z{ ~ r(t)x2) ■ 


Q 2 


V A{t) 


In the above nonisospectral eigenvalue problem, spectral parameter £( t ) obeys the following equation 

C(t) = 

where /i is a hidden complex constant and T(t) is an arbitrary function of time and 

r (*) = 


a{t) = 


A{ty 


a {tf = r(t ) 2 - r \t) 


(24) 

(25) 

(26) 

(27) 

(28) 


It is known [22[ [23] that the coupled GP equation is completely integrable only if the transient trap \{t) and the 
scattering length a(t ) are connected by the following equation (obtained by substituting eqn.(26) and eqn.(|27|) into 
eqn.(28)) 


1 d 2 a(t) + 1 fda(t)\ 2 


2 a(t) dt 2 a 2 (t ) 




(29) 


To generate dark vector solitons of the above GP eq. ([20]) , we choose the following non zero plane wave solution as 
the seed 


q. — c . e [i(0'iX-(al /2+<ncf)t)] 


( 30 ) 
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FIG. 1: The density profiles of the dark solitons in the presence of periodic modulated potential for the choice of parameters 
T(t) — cos(ujt + S)e at , a;=6,<5=0,<7=0.4,=0.3, <aio=0.1, /3io=0.3, xi —0.1, (5i=0.2,ai=-a2=l,ci=C2=l 



FIG. 2: Contour plots of fig [I] exhibiting the beating effect in the propagation dark solitons 



FIG. 3: Collisional dynamics of dark solitons for the choice of parameters a(t) = 0.1, T(t) = 0.1 x 10 2£, a\ — 0.1, c ^2 = 0.25, 
Pi = 0.3, p 2 = 0.25, <5i = 0.1, <5 2 = 0.2, X i = 0.3, X2 = 0.4, = 0.85*, e{ 2) = 0.5 such that | 4 j) | 2 + \e ( ^\ 2 = 1, (j = 1, 2),oi = 

— Ci2 = 1, Cl = C2 = 1 































where 


4 ? 

4 ? 

4 ? 


_ e -2iC(t)s-6iC(t) 2 t^ 
_ e i((t) x +3i((t) 2 t ^ 

_ p i((t)x+3i((t) 2 t 


Employing the gauge transformation method, one obtains the dark soliton solution for the coupled GP equation (20) 
(for the defocussing case keeping the temporal scattering length a(t)=-a(t)) of the following form 


V’h = \j ^y2c i e^ 1) ^ 1 ( t)tanh(6i)e l{ - t 1+ fW x 2 \ 

\j ^t) 2 c * 4 1) ^ 1 (t)tanh(0i)e l ^ 1+fWs ^\ 


4 1] 


where 


and 


0i = 2/3ix + 4 J ai/3idt — 25i, 

= 2aix + 2 J ( a i ~ 0i)dt + A — 2%i, 

A = e [^(aiiE-(aJ/2+o-iC?)t)] 


(32) 

(33) 

(34) 

(35) 

(36) 


with aq = ceioe^ r ^ dt l, /?i = j3 10 e^ r ^ dt ^ while 8\ and xi are arbitrary parameters and e^\e^ are coupling 
parameters connected by the relation | 2 + l^^l 2 = 1? U = 1?2). The dark solitons given by eqs.(32) and |33j ) are 
identical to the one reported in [23 . The density profile of dark solitons given by eqs.(32) and (33) is shown in fig. 0 
while its contour plots shown in fig. 0 display the beating effect of the dark solitons during time evolution. It should 
be mentioned that this beating effect arises due to the temporal nature of the harmonic trap |23 , [24] . Eventhough the 
nature of the solitons depends on the scattering length a(t) and trap frequency A (t) (or T(t)) in accordance with the 
integrabilily condition given by equation. (29), the trap frequency T(t) predominates over the scattering length. As we 
have chosen the trap frequency T(t) to be a periodic wave with exponentially varying amplitude (a > 0), the amplitude 
of the soliton visibly gets higher for large ”t” as it is evident from figs jl] and fig. (J2|. This beating effect of dark solitons 
is consistent with the experimental results discussed in ref. [25]. Our paper also gives a simple experimental protocol to 
observe ” beating effect” in the collisional dynamics of dark solitons without the impact of external thermal cloud.The 
collisional dynamics of two solitons is shown in fig. One can easily extend this modified gauge transformation 
approach to construct multi dark soliton solutions and study their collisional dynamics. 

It should also be mentioned that the present integrable model does not have the luxury of observing parametric 
resonance excitations in dark solitons (BECs) by virtue of the constraint imposed by the integrability condition given 
by eq.(|29|) as one may not be able to choose the frequency of the trap T(t) and a(t) desirably. 


(ii) Coupled Nonlinear Schrodinger (CNLS) equation 


Under the dependent variable transformation 


Vh = 


1 


-^> 1 (X,T)exp(i[T(t)x 2 }) 


^ = y^W $2(x ’ T)e ^ [r( ^ 2]) 


(37) 

(38) 


where is an arbitrary function with spatial and temporal variables chosen as X = and T = l 2 with T(t), A (t) 
given by 


r t + 4r 2 + a 2 = o j t - 4r i = o, 


(39) 
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1 


FIG. 4: Density profiles of dark solitons of coupled NLS equation for the choice — 0.3, a±o = 0.3, /3\o = 0.5, Xi — 0.3, Si = 
0.5,ai — —a 2 — l,c\ — C 2 — \ 



FIG. 5: contour plots of fig. 0 showing the trajectory of dark solitons 




eqn.(20) can be reduced to the standard coupled nonlinear Schodinger equation in the normal GVD region of the 
following form 


i<&i T — <&ixx + Kl^il 2 + 1^2 | 2 )^i — 0, (40) 

i<&2T — &2XX + ^(l^il 2 + |^2 | 2 ) < ^ > 2 = 0. (41) 


The above coupled nonlinear equation is nothing but the celebrated integrable model proposed by Manakov m- 
The dark soliton solution of the CNLS equations employing the gauge transformation approach is of the following 
form (defocussing- defocussing Manakov model for v = — 1) 


where 


and 


= —ci£^/3itanh(6i)e' l ( 

= —C2£^Pitanh(6i)e L ^-^ 1 \ 

6i = 2p\x + 8oq/3i£ — 2S\, 

fi = 2ol\x + 4(af — j3l)t + A — 2%i, 


\ = e [i(aiX-(a?/2+aic?)t)] 


(42a) 

(42b) 


(43) 


subject to \£i\ 2 + \ e 2 ^\ 2 = F U = 1?2) where represent coupling parameters, otj = Re[/j] and /3j = Im[p\ while 
j — 1,2 are the real constants and fi is the hidden spectral parameter. The dark soliton given by eqn.([42| is identical 
to the one reported in ref [26]. The density profile of dark solitons and its trajectory are shown in figs. 0 and 0 
respectively. The collisional dynamics of dark solitons of the Manakov model is displayed in fig. 0. 
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FIG. 6 : Collisional dynamics of dark solitons of the coupled NLS equation for the parametric choice ai = 0.1, a 2 = 0.25, 
Pi = 0.3, p 2 = 0 . 2 , ft = 0 . 1 , 62 = 0 . 2 , xi = 0.3, X 2 = 0.4, e™ = 0.85i, ef 0 = 0.5 such that \s[ j) \ 2 + \e ( p\ 2 = l/v,(j = 1 , 2 ), 

Cll — —Cl2 = 1, Cl = C2 = 1 


IV. CONCLUSION 

In this paper, we have formulated a simple algebraic approach by harnessing the freedom in the celebrated gauge 
transformation approach to construct dark solitons of coupled NLS type equations. As an application, we have 
constructed the dark solitons of coupled GP and coupled NLS equations and studied their properties. It should be 
emphasized that the present approach is purely algebraic and enables one to generate multi dark soliton solutions from 
a trivial/nontrivial seed solution and it is quite superior to other analytical techniques like Darboux transformation 
method [27] . Hirota method [26[ [28] and 1ST [29]. As far the limitations of this approach are concerned, it should 
be emphasized that this method generates only bright and dark solitons for zero and non zero seeds respectively. We 
have not yet exploited it to generate breathers, Ma solitons, rogue waves amongst others. 
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